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ABSTRACT 

We present an anisotropic analysis of the baryonic acoustic oscillation (BAO) scale in the 
twelfth and final data release of the Baryonic Oscillation Spectroscopic Survey (BOSS). We 
independently analyse the LOWZ and CMASS galaxy samples: the LOWZ sample contains 
contains 361 762 galaxies with an effective redshift of zlowz = 0.32; the CMASS sample 
consists of 777 202 galaxies with an effective redshift of zcmass = 0.57. We extract the BAO 
peak position from the monopole power spectrum moment, ao, and from the p? moment, a 2 , 
where pi is the cosine of the angle to the line-of-sight. The /x^-moment provides equivalent 
information to that available in the quadrupole but is simpler to analyse. After applying a re¬ 
construction algorithm to reduce the BAO suppression by bulk motions, we measure the BAO 
peak position in the monopole and ^^-moment, which are related to radial and angular shifts in 

scale. We report i?(zLowz)?'s( 2 d) = (11.60± 0.60) • 10^ kms“^ and i:)A( 2 ^LOWz)/?’s( 2 d) = 

6.66 ± 0.16 with a cross-correlation coefficient of thda — 0.41, for the LOWZ sample; and 
H^zcMASsysizd) = (14.56 ± 0.37) • 10^ kms"! and Da{zcmass)/ rs{zd) = 9.42 ± 0.13 
with a cross-correlation coefficient of thda = 0.47, for the CMASS sample. We demon¬ 
strate that our results are not affected by the fiducial cosmology assumed for the analysis. 
We combine these results with the measurements of the BAO peak position in the monopole 
and quadrupole correlation function of the same dataset ( [Cuesta et al.|2016| companion pa¬ 
per) and report the consensus values; i7(zLOWz)?'s(.Zd) = (11.63 ± 0.69) • 10^kms“^ 
and DA{zi^o^z)lfa{zd) = 6.67 ± 0.15 with thda = 0-35 for the LOWZ sample; 
H{zcMASs)ra{zd) = (14.67 ± 0.42) • lO^kms”! and DA{zcyiASs)/rs(zd) = 9.47 ± 0.12 
with thda = 0-52 for the CMASS sample. 

Key words: cosmology; theory - cosmology; cosmological parameters - cosmology; large- 
scale structure of Universe - galaxies; haloes 
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1 INTRODUCTION 


2 DATA AND MOCKS 


The Baryon Acoustic Oscillation (BAO) signal in the clustering 
of galaxies provides a robust route to measure the cosmological 
expansion rate ( [Eisenstein et al.|2005| >. The angular separation of 
galaxies conveys different cosmological information compared to 
the radial separation. The observed projection of the BAO-scale 
in the angular direction depends on the angular diameter distance, 
while the radial projection is a function of the Hubble expan- 
sior0 The anisotropic information can be extracted from measur¬ 
ing the BAO peak position in moments of the correlation function 
or Fourier-space power spectrum. 

The Baryon Oscillation Spectroscopic Survey (BOSS; pZ)aw-| 
|son et al.|2013| ), which is part of SDSS-III ( [Eisenstein et al.|201 1 1 , 
has provided the largest set of numbers of spectroscopic galaxy ob¬ 
servations made to date; in this paper we analyse the final. Data 
Release 12 (DR12; [Alam et ar]|2015| ( sample. Our work follows 
on from measurements made for DR 11 samples presented in |An-| 
jderson et al.j ( [2014^ and [Beutler et al.j j2014^ . We analyse the 
low-redshift (LOWZ) and high-redshift (CMASS) samples inde- 
pendently, which were targeted using different algorithms ([Reidj 


2.1 The SDSS III BOSS data 


et al.||20f5 i. DR12 contains approximately 1000 deg^ more solid 


angle than DRl 1, an increase of approximately 12%. Additionally, 
the methodology adopted to create galaxy catalogues has been im¬ 
proved jReid et al.|2015[[Ross, Percival & Manera|2015| l, leading 
to improved understanding, and mitigation, of potential systematic 
errors. 

In this paper we present a Fourier-space analysis of the line-of- 
sight (LOS) dependent clustering, extracting the BAO position. The 
monopole and quadrupole calculated were analysed to measure the 
Redshift-Space Distortion (RSD) signal in jGil-Marm et al.| ( |2015| l, 
hereafter Paper I. In this paper we present a new method for fitting 
the -moment of the power spectrum rather than the quadrupole 
to obtain the anisotropic information. As the /i^-moment is a sim¬ 
ple linear combination of monopole and quadrupole, this change 
is lossless. Fitting to the ^^-moment allows the use of the same 
equations as for the monopole and quadrupole (with different pa¬ 
rameters) simplifying the modelling, and particularly the way in 
which we can isolate BAO information from the broadband sig¬ 
nal. Note that using the ^^-moment instead of the quadrupole is a 
mere fitting technique that simplifies the identification of the BAO 
peak in the higher order multipoles of the power spectrum. How¬ 
ever, one can still use the estimators used in Paper I to compute the 
monopole and quadrupole, and reconvert them into monopole and 
^^-moment. Since the /r^-moment is just a linear combination of 
the monopole and quadrupole, the result is identically the same to 
measure the -moment using the estimator described in 
|eni| ( |20T5) . 

The outline of our paper is as follows: In we introduce the 
data, the mocks used to compute the covariance matrices and test 
for systematics, and briefly describe the reconstruction technique. 
In we present the formalism used to calculate the monopole 
and /i^-moments. The method applied to fit the BAO is described 
in Section]^ The measurements of the BAO peak position in the 
DR12 BOSS galaxies are presented in for both pre- and post- 
reconstructed catalogues. In i|^we test for potential systematics of 
the model and of the fitting process using the post-reconstructed 
mocks. Finally ^presents the conclusions of this paper. 


Bianchi 


Smee et al. 


As part of the Sloan Digital Sky Survey III ( [Eisenstein et al.|201 1 1 
the Baryon Oscillations Spectroscopic Survey (BOSS) ( Dawson| 
e^^alJ20^ I measured spectroscopic redshifts jBolton et al.|2012 


2013^ for more than 1 million galaxies and over 200 000 


quasars. The galaxies were selected from multi-colour SDSS imag¬ 


ing ([Fukugita et al.[[1996| [Gunn et al.[[1998[ [Smith et al.[[2002[ 
[Gunn et al.|2006[ [Doi et al.|2010t focussing on the redshift range 
of 0.15 ^ 0.70. The galaxy survey used two primary target 

algorithms, selecting samples called LOWZ, with 361762 galax¬ 
ies in the final data release DR 12 ( [Alam et al.[[2015) between 
0.15 sC 2 s: 0.43 and CMASS, with 777 202 galaxies in DR12 
between 0.43 ^ z Si 0.70. The full targeting algorithms used and 
the method for calculating the galaxy and random catalogues are 
presented in [Reid et al.[p015) l. The samples jointly cover a large 
cosmic volume (Uefi = 7.4 Gpc^) with a number density of galax¬ 
ies that ensures that the shot noise does not dominate at BAO scales 
at the relevant redshifts. Obviously in the edges of the redshift-bin 
of the LOWZ and CMASS samples the shot noise has a more rel¬ 
evant role than close to the center of the redshift-bin, just because 
of the difference in the number density of galaxies. However, the 
power spectrum signal is dominated by those galaxies of the center 
of the bin, which is less affected by shot noise. Full details of the 
catalogues are provided in [Reid et al.[j2015) . 

In order to correct for several observational artefacts in the 
catalogues, the CMASS and LOWZ samples incorporate a set of 
weights: a redshift failure weight, Wy.{, a fibre collision weight, trifc, 
and a systematic weight, Wsys (CMASS only), which combines a 
seeing condition weight and a stellar weight ( [Ross et al.[2012[[An-[ 
[derson et al.[[2014^ [Reid et ^[2015^ . Hence, each galaxy target 
contributes to our estimate of the target galaxy density field by 


Wc = msys(mrf + Wic — 1). 


( 1 ) 


The redshift failure weights account for galaxies that have been ob¬ 
served, but whose redshifts have not been measured: nearby galax¬ 
ies, which are approximated as being “equivalent”, are up-weighted 
to remove any bias in the resulting field. The fibre collision weight 
similarly corrects for galaxies that could not be observed as there 
was another target within 62", a physical limitation of the spectro¬ 
graph (see [Ross et al.[j2012^ for details). The systematic weight ac¬ 
counts for fluctuations in the target density caused by changes in the 
observational efficiency. This effect is only present for the CMASS 
sample, which relies on deeper imaging data; such a weight is not 
required for the brighter LOWZ sample ( [Tojeiro et al.[2014) . 

Additionally, we adopt the standard weight to balance regions 
of high and low density ([Feldman, Kaiser & Peacock[1994|[Beutler[ 
[et al.[2014[ l. 


wpKp(r) = 


msys(r) 


( 2 ) 


Wsya(r) -I- Wc(r)n(r)Po ’ 

where n is the mean number density of galaxies and Pq is the am¬ 
plitude of the galaxy power spectrum at the scale where the er¬ 
ror is minimised. We assume Po = 10 000 h~^ Mpc®, which cor¬ 
responds to the amplitude of the power spectrum at scales k ~ 


0.10/iMpc I Reid et al. 2015 I. 


^ Negligible evolution is expected over the redshift interval between galax¬ 
ies separated by the BAO scale 


2.2 The mock survey catalogues 

Mock samples are a key component in the analysis of precision 
cosmological data provided by galaxy surveys. They are a fun- 
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Da 

Clbh? 

h 

Tls 

Fid. 

0.31 

0.69 

0.023569 

0.7 

0.9624 

MD-Patchy 

0.307115 

0.692885 

0.022045 

0 . 6 m 

0.96 

QPM 

0.29 

0.71 

0.022442 

0.7 

0.97 


Table 1. Cosmological parameters chosen for the fiducial cosmology, for 
the MD-Patchy cosmology and QPM cosmology. 


damental requirement to test and analyse the large-scale structure 
and they help determine systematic errors on the measurements. 
Most of the relevant large-scale physics is captured by approxi¬ 
mate methods, so we do not necessarily need to base mock cat¬ 
alogues on full N-body cosmological simulations; small numbers 
of N-body simulations can instead be used to calibrate a more ef¬ 
ficient scheme. In this paper we use mocks created by two dif¬ 
ferent approaches: MultiDark-Patchy BOSS DR12 mock^ (here¬ 
after MD-Patchy mocks) ( [Kitaura et al. 2016 companion paper) 
and Quick-Particle-Mesh mocks |White, Tinker & McBride|2014t , 
hereafter QPM mocks. Both schemes incorporate observational ef¬ 
fects including the survey selection window, veto mask and fiber 
collisions. In Paper I we demonstrated that the covariance matrices 
obtained using these two types of mocks were similar and yielded 
similar errors-bars and cross covariance parameters (see Fig. 3 and 
13 in Paper I). Because of this, in this paper we perform two parallel 
analyses of the data using MD-Patchy and QPM mocks. Since the 
measurements and the errors coming from the two sets of mocks 
are very similar, we average both the measurements and errors in 
order to obtain a single value for each parameter. Since the covari¬ 
ance matrix is estimated from a set of mocks, its inverse is biased 
due to the limited number of realizations. We account for this effect 
by applying the correction proposed by |Hartlap, Simon & Schn^ 
|der| j2007| l. In addition to this scaling, we have to propagate the 
error in the covariance matrix to the error on the estimated param¬ 
eters. We can do this by scaling the variance for each parameter by 
the factor of eq. 18 of |Percival et al.| ( f2014^ . The corrections to the 
error bars described above are small (< 2% for QPM and < 1% for 
MD-Patchy), but are included in all the error bars quoted in this 
paper. 


2.3 Fiducial Cosmology 

In this paper we analyse the data assuming a fiducial cosmological 
model. The values of this cosmological model, as well as the cos¬ 
mologies of QPM and MD-Patchy mocks, are presented in Table 

□ 

In ^we show that the arbitrary choice of cosmology has a 
systematic effect of < 0.3% on the peak position values obtained 
from the QPM mocks. We have checked the effect of this systematic 
error by adding it in quadrature along with the statistical errors for 
the data. We have found that this represents just a 5% increase of 
the size of the total error-bars budget for the case of the position 
of the BAO peak in the isotropic post-recon signal in the CMASS 
sample, and much less in the rest of the variables. Because of this 
we consider this effect as sub-dominant and we will not consider it 
in the results of this paper. We present a further discussion of this 
point later in 


^ http://data.sdss3.org/datamodel/index-files.html 
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2.4 Reconstruction 

The initial BAO signal in the clustering is damped by the comov¬ 
ing motions of galaxies, potentially reducing the fidelity of BAO 
measurements. Using the observed density field, we can model the 
galaxy motions, and move the over-density back to its “original” 
position, recovering a fraction of this signal ( [Eisenstein et al.|2007| ). 
The fraction of the signal recovered is dependent on the density of 
galaxies, which limits how well the displacement field can be de¬ 
termined ( [Burden et al.||2014t . Redshift-Space Distortions (RSD) 
complicate the implementation of this method, as they produce a 
LOS-dependent distortion that has to be included when estimating 
the displacement field. In Fourier-space, an iterative method can 
be used to compensate for this effect ([Burden, Percival & Howlett| 
[20T5| |, where a correction is computed at each step using the prior 
estimate of displacements. An alternative approach is to grid the 
galaxy over-density field, and determine RSD and displacements 
using a finite-difference routine based on the values of the over¬ 
densities at grid-points ( jPadmanabhan, White & Cohn|2009| l. Both 
methods provide consistent results; we use the finite-difference 
method here. We will make measurements both before applying 
this algorithm (pre-recon), and after (post-recon). 


3 POWER SPECTRUM ESTIMATOR 

We build the galaxy power spectrum multipole estimator by defin¬ 
ing the function 

Ff) = [wc(r)w(r) - a^ann^r)], (3) 

where n and Ua are, respectively, the observed number density of 
galaxies and the number density of a synthetic catalog Poisson sam¬ 
pled with the same mask and selection function as the survey with 
no other cosmological correlations. The functions Wc and wfkp 
were defined in Eqs. (jTJ and ^ respectively. The factor Qran is 
the ratio between the weighted number of observed galaxies over 
the random catalogue, aran = Wc/Na, where Na denotes 

the number of objects in the synthetic catalog and Agai the number 
of galaxies in the real catalog. The factor I 2 normalises the ampli¬ 
tude of the observed power in accordance with its definition in a 
universe with no survey selection. 

We compute the power spectrum multipoles using the imple¬ 
mentation of the Yamamoto estimator ^Yamamoto et al.|2006 |l pre¬ 
sented in jBianchi et al.| ( |2015| l, which is based on using multiple 
Fast Fourier Transforms, keeping the relevant LOS information by 
approximating the LOS of each pair of galaxies with the LOS of 
one of the two galaxies (see section 3 of Paper I for more details on 
the algorithm used). This is a reliable approximation on the scales 
of interest, which clearly improves on assuming a single fixed LOS 
for the whole survey for the quadrupole, but will eventually break 
down at large scales and high order multipoles ( |Yoo & Seljak | 2015} 
Samushia, Branchini & Percival|2015| l. 

We use a random catalogue of number density of na{v) = 
c*rand(r), with — 50, for pre-recon, and — 20 for post- 
recon catalogues. We place the LOWZ and CMASS galaxy samples 
on 1024^ quadrangular grids, of a box of side Lh = 2300 /i“^Mpc 
for the LOWZ galaxies, and Lb = 3500 /i“'^Mpc to fit the CMASS 
galaxies. This approach corresponds to a grid-cell resolution of 
3.42 /i“^Mpc for the CMASS galaxies and 2.25 /i“^Mpc for the 
LOWZ galaxies. The fundamental wave-lengths are therefore kj = 
1.795 • 10“® /iMpc"^ and kf = 2.732 • 10“^ /iMpc"^ for the 
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CMASS and LOWZ galaxies, respectively. We apply the Cloud-in- 
Cells scheme (CiC) to associate galaxies to grid-cells, and bin the 
power spectrum fe—modes in 60 bins between the fundamental fre¬ 
quency kf and a maximum frequency of = 0.33 /i Mpc, with 

width Alogio * = [logio(fcM) - logio(fc/)] / 60 . 

We limit the scales fitted as follows: our procedure for deter¬ 
mining the largest scale for the fitting process is based on limit¬ 
ing the impact of the systematic weights, and is presented in ap¬ 
pendix A of Paper I. We limit scales to k > 0.02for 
the monopole and k > 0.04 Mpc“^ for the quadrupole. The 
smaller scale used for the BAO peak determination is fcmax = 
0.3 fe Mpc“^, as for smaller scales the BAO information is quite 
limited. 

In this paper we work with the combined north and south 
samples, NGC-l-SGC, for all the power spectrum multipoles. This 
combination is performed by averaging both NGC and SGC com¬ 
ponents weighted by their area, Pngc+sgc = (Pngg^ngc + 
PsgcA-sgc)/ (j4ngc + j4sgc) 5 both for LOWZ and CMASS sam¬ 
ples. The values of the areas are: A^g^^ = 5836 deg^, = 

2501 deg^ = 6851 deg^ and = 2525 deg^ 


4 FITTING THE POWER SPECTRUM MOMENTS 

4.1 Modelling the power spectrum moments 

We model the power spectrum multipoles in order to measure 
the position of the BAO features and marginalise over the broad¬ 
band information. To reduce the complexity of the fits, we choose 
to fit the monopole and the ^^-moment, defined as ' = 
2/5P^^^ (k) + P^°^ (k), rather than the monopole and quadrupole. 
As the pair monopole - -moment are a linear transformation of 
the pair monopole - quadrupole, the information content is the same 
in both. P^'^ ^ has the property that, in linear theory, the shape is 
the same as P^'^\ and we can therefore use the same modelling pro¬ 
cedure for both, albeit with different free parameters describing the 
deviations from linear theory to account for LOS dependent effects. 

To create the model moments, we start by computing the lin¬ 
ear power spectrum at a given cosmology, Pun, which we gener¬ 
ate using Game ( [Lewis, Challinor & Lasenb^|2000^ . The linear 
power spectrum was separated into two components, one contain¬ 
ing the BAO oscillations, Oun, and a smooth component, Piin.am; 
Piin(fc) = Oiin(fe)Piin,sm(fe). TMs Separation was performed us¬ 
ing the same method applied to the data, but employing an analytic 
model for the BAO and smooth model based on the fitting formu¬ 
lae of lEisenstein & Hu| ( |1998) . The position of the BAO peak is 
described by a, which parametrises the features of the oscillations 
as a function of k in the Oun-function, Oiin{k/a). We use a su¬ 
perscript “(0)” for the monopole and “(2)” for the value measured 
from the -moment. We adopt the model described in 
|et al.| ( [20T4l l to account for deviations in the smooth fit away from 
this linear model. The full model is 

f’model(fe;a) = f’model,sin(fe) |l + [Olin{k / a) - | , 

(4) 

where Pmodei.sm is a phenomenological parametrisation of the 
non-linear power spectrum monopole with no-BAO, 

PmodeMm(fc) =B"Plin,sm(A:)+Alfc-fA2 + ^ + ^ + ^, (5) 

and B, Ai and Eni are free parameters that account for redshift 
space distortions and nonlinearities. Together with a, there are 8 
free parameters for each model. These parameters are allowed to 


Anderson 


be different in the fits to the monopole and moment and we 
adopt superscripts “(0)” for the monopole and “(2)” for the fit to 
the /r^-moment, where appropriate. 

The survey window acts as a convolution of the power spec¬ 
trum moments, smoothing the true power spectrum to produce that 
observed. As the survey window affects the BAO, it cannot be de¬ 
scribed by a variation of the free parameters described above. To in¬ 
clude the window effects, we follow the approach described in §5.4 
of Paper I, producing matrices to account for these effects through 
a discrete convolution 

PitXki) = 

j 3 

PH.(ki) = 

3 3 

(6) 

where, Wj"" are the elements of the window-survey matrix, 
P(") model are the pre-masked models described by Eq. (|^, and 
^w'in products for the observed moments. Details of 

how the Wj"" elements are estimated can be found in eq. (14) of 
Paper I. This model is able to fit both pre- and post-recon power 
spectrum moments with residuals that are of far lower significance 
than the BAO signal (see ©• 


4.2 Finding the best-fitting parameters 

The monopole and -moment are assumed to be drawn from a 
multi-variate Gaussian distribution, whose covariance can be com¬ 
puted using the galaxy mocks described in §2.2| For the pre-recon 
analysis we consider two different sets of mocks, QPM and MD- 
Patchy, whereas for the post-recon analysis we only use QPM 
mocks. For details about the computation of the covariance ma¬ 
trices we refer the reader to section 6.1 of Paper I. 

Eq.(|4} includes 8 free parameters for each of the two power 
spectrum moments, so we have 16 free parameters in total to be fit¬ 
ted to each set of measurements. To investigate this parameter space 
we use a bespoke Markov-Chain Monte-Carlo (MCMC) routine, 
run for each fit for 10® steps split into 10 sub-chains, which satisfies 
the convergence check that each sub-chain gives consistent results. 
Each chain was started close to the best-fitting locations found us¬ 
ing a downhill simplex routine, and a bum-in of 10® steps was re¬ 
moved from each. To reduce the computational burden, each chain 
was thinned by a factor 10^ before measuring marginalised param¬ 
eters and errors. We only consider the interval 0.8 < Qi < 1.2 
for Qo and a2 to avoid unphysical solutions and to reduce the con¬ 
vergence time of the MCMC routine. This is a very conservative 
prior which is never hit by the MCMC routine, and that happens 
to be at more than 5a of the final results. Due to noise, secondary 
maxima where the BAO signal is located within the large-scale, 
noisy, part of the power spectrum, can cause a mobility issue for 
the chains. Given that the data storage and analysis are computa¬ 
tionally expensive, we choose to run relatively long chains, which 
are less affected by the secondary maxima, and thin them to allow 
the subsequent analysis to be performed quickly. 

The correlation among the parameters of interest, ao and 0:2 is 
shown later in Fig.[^ The correlation among the a-parameters and 
the rest of nuisance parameters, B, Ai and Eni, is found to be very 
weak or consistent with 0 (|r| < 0.10). The strongest is a weak 
correlation of r ~ 0.27 between ao and Eni; and between a.2 and 
Eni for the CMASS sample; and r ~ —0.30 between ao and Eni; 
and between 02 and Eni for the LOWZ sample. 
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4.3 Interpreting the measured BAO scales 


Our BAO scale measurements are given by ao and 02 , and their 
covariance. To translate these parameters into easily-to-model cos¬ 
mological measurements, we adopt the standard definitions of the 
Alcock-Paczynski (AP) parameters, 


^ ^ DA{z)r^s\za) 

H{z)rs{zd) ’ D^{z)rs{zd) ’ 


(7) 


Our measurements of ao and 0:2 each provide degenerate measure¬ 
ment of a|| and a^. Ross, Percival & Manera ( 2015|( showed that 
we should expect a degeneracy of the kind. 


m + n 


m n 

ail ax, 


( 8 ) 


where for the post-recon monopole (and RSD removal), m = 1/3 
and n = 2/3, and for the /r^-moment, m = 3/5 and n = ‘Ij'o. 
Our window function convolution accounts for any deviations from 
these ideal solutions caused by the survey geometry. Thus, the AP 
parameters, ay and ax can be written in terms of the position of 
the BAO peak in the monopole and -moment. 


- 3/2 5/2 
— a, ' n, ! ■ 


0 


a. 


ax = a, 




(9) 


Using these expressions we relate the BAO peak position in the 
monopole and -moment to the Hubble parameter and the angular 
distance as, 

H{z)rs{zd) = [H{z)rs{zd)f’^af‘^a'^^^'^ ( 10 ) 

and 

DA{z)/rs{zd) = [DA{z)/rs{zd)f'^ay'^Q.2^^'^. ( 11 ) 

For the fiducial cosmological model of Table 
[JT(2Lowz)fs(2d)]‘^‘^ = 11.914 • 10®kms“^ and 

[i9A(2LOWz)/fs(2d)]*’‘^ = 6.667 for the LOWZ sam¬ 

ple and [JT(2cMASs)fs(2d)]*’‘^ = 13.827 ■ lO^kms”^ and 
[i9A(2CMASs)/t'3(2d)]®'^ = 9.330 for the LOWZ sample. In this 
fiducial cosmological model we have rs{zd) = 143.70Mpc. We 
report these parameters and its correlation coefficient as a final 
result. 


5 TESTS ON GALAXY MOCKS 

In this section we test for potential systematics in the BAO model as 
well as those associated to the fitting algorithm. We focus on QPM 
post-recon mocks, which we analyse assuming using two different 
cosmologies (with Q.m = 0.29 or 0^ = 0.31) and two different 
BAO models (with Q.bh^ = 0.022442 or = 0.023569). For 
the four different combinations of 0^ and flbh^, we test that we 
are able to recover the expected BAO shifts positions even when 
the cosmology or the BAO model assumed for analysing the mocks 
does not necessarily match their true values. In the case where both 
the cosmology and the BAO model assumed match those values 
from the mocks, the expected values for ao and 02 are 1. However, 
when either the cosmology or the BAO model assumed differ from 
those of the mocks, we expect that the values for ao and a 2 will 
differ from 1, with the expected values given by Eqs. |^. We will 
refer to these expected BAO peak position as and 

We start by computing the best-fitting values of ao and 02 for 
the average monopole and /i^-moment over the 1000 realisations of 
the mocks, for the four different combination of Qm and Qbtd (la¬ 
belled a{P) in Table|^. By averaging the power spectrum moment 
over the realisations we reduce the noise-dependent systematic er¬ 
rors to a negligible level compared to the systematic errors of the 


model. Therefore, when we compare the obtained a; with its ex¬ 
pected value we are exclusively testing the potential systematics of 
the BAO model itself. Given the results of Tablel^associated to the 
statistic a{P) we can state that the BAO model is able to reproduce 
the expected value of the BAO peak position in both the monopole 
and /r^-moment with < 0.30% accuracy. 

In order to test for noise-dependent systematic errors we have 
also computed expected value for ao and a 2 for each mock and 
average among them. Unlike the former case, now we do not can¬ 
cel the noise-dependent systematic associated to each individual fit. 
The results are listed in Table|^ labelled a{P) (i.e. the average of 
the BAO peak position of the mocks). For some values of the cos¬ 
mology and BAO model we observe that the deviation between the 
expected and the recovered values has significantly changed with 
respect to fitting the average power from the mocks, a{P). This 
result is due to the remaining statistical noise associated to each in¬ 
dividual mock fit. In general, we are able to recover the BAO peak 
position in both the monopole and ^^-moment with ^ 0.25% ac¬ 
curacy. These values are significantly below the statistical errors on 
the measurements, in a similar way as we have seen for the a{P) 
statistic. We therefore conclude that the potential systematic errors, 
both noise-dependent and model-dependent, do not play any im¬ 
portant role in the BAO peak position estimation in the monopole 
or /r^-moment, where the statistical errors of the survey dominate. 
We therefore do not apply any correction to the obtained results 
from the data. 


6 RESULTS 

Fig. □ presents the post-recon power spectrum monopole (blue 
squares), quadrupole (blue circles) and /r^-moment (green trian¬ 
gles), for the LOWZ and CMASS samples of the DR12 data mea¬ 
surements (top and bottom panels as labeled) from the combi¬ 
nation of the NGC and SCG. The coloured lines show the best¬ 
fitting model given by Eq. with the BAO peak position for the 
monopole (ao) and for the /r -moment ( 02 ). For simplicity we only 
show the best-fitting model according to the covariance matrix ex¬ 
tracted from the QPM mocks. We will later show in Table the 
results from the data assuming both QPM and MD-Patchy co- 
variance. The model for the quadrupole has been obtained from 
the models for the monopole and -moment. As the BAO peak 
positions in the monopole and -moment are similar to those 
in the fiducial model, there is no apparent BAO signature in the 
quadrupole. Within each panel, we also display the data and best¬ 
fitting model for the monopole and -moment divided by the 
smoothed power spectrum P^m of Eq. § for the best-fitting model. 
For clarity, the bottom sub-panels show the residuals between the 
measurement and the model, AP = _ pdata^ divided by 

the la error of the data. The black dashed lines marks the Icr and 
2a deviations. The model is able to accurately reproduce the BAO 
features both in the monopole and -moment for both LOWZ 
and CMASS sample. As expected, the BAO features are more sig¬ 
nificant in the monopole than in the -moment, as the former 
has a higher signal-to-noise ratio. Similarly, the BAO features are 
stronger in the CMASS sample than in the LOWZ. 

At large scales, the quadrupole measured for both LOWZ and 
CMASS samples is small, consistent with the reconstruction pro¬ 
cess removing the linear component of the RSD. However, the re¬ 
construction process is not able to fully suppress the whole RSD 
signal and non-linear components are left, as is the AP effect caused 
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Figure 1. The measured LOWZ (top panel) and CMASS (bottom panel) DR12 post-recon, monopole (blue squares), quadrupole (red circles) and /i^-moment 
(green triangles) power spectra. For all the cases the measurements correspond to a combination of the northern and southern galaxy caps according to 
their effective areas as described in The error-bars are calculated from the dispersion of measurements using the QPM mocks. The red, blue and green 
lines correspond to the best-fitting model of Eq. with the BAO peak position as a free parameter. Within each panel we also present the power spectrum 
monopole and /i.^-moment divided by the smooth power spectrum calculated in our fit to the data. For the monopole and -moment we see how the model 
is able to capture the BAO features observed in the data. For clarity, the bottom sub-panels show the residuals between the measurement and the model, 
AP = p^odel _ pdata^ divided by the Icr error of the data. The black dashed lines marks the la and 2a deviations. 
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BAO measurement from the LOS-dependent power spectrum ofDRll BOSS galaxies 


1 


Statistic 

Sample 

BAO model 

Cosmology 

exp 

ao - «o 

“0 

exp 

OL2-a2 

ar'’ 


LOWZ 


= 0.023569 


= 0.29 

= 0.31 

0.00137 ± 0.00074 
0.00061 ± 0.00048 

0.976687 

0.982765 

0.0022 ± 0.0010 

0.00103 ± 0.00066 

0.976687 

0.984004 

a{P) ■ 


= 0.022442 


= 0.29 

= 0.31 

0.00137 ± 0.00048 
0.00131 ± 0.00050 

1.000000 

1.006223 

0.00175 ± 0.00068 
0.00158 ± 0.00069 

1.000000 

1.007491 

CMASS 


= 0.023569 


= 0.29 

= 0.31 

-0.00253 ± 0.00031 
-0.00274 ± 0.00031 

0.976687 

0.986808 

-0.00257 ± 0.00044 
-0.00290 ± 0.00044 

0.976687 

0.988829 



= 0.022442 


= 0.29 

= 0.31 

-0.00098 ± 0.00032 
-0.00127 ±0.00032 

1.000000 

1.010363 

-0.00102 ± 0.00044 
-0.00147 ± 0.00045 

1.000000 

1.012431 


LOWZ 


= 0.023569 


= 0.29 

= 0.31 

0.00173 ± 0.00050 
0.00204 ± 0.00051 

0.976687 

0.982765 

0.00227 ± 0.00071 
0.00254 ± 0.00073 

0.976687 

0.984004 

5(P) . 


= 0.022442 


= 0.29 

= 0.31 

0.00153 ±0.00051 
0.00177 ± 0.00052 

1.000000 

1.006223 

0.00140 ± 0.00072 
0.00155 ± 0.00073 

1.000000 

1.007491 

CMASS 


= 0.023569 


= 0.29 

= 0.31 

-0.00243 ± 0.00032 
-0.00224 ± 0.00033 

0.976687 

0.986808 

-0.00244 ± 0.00045 
-0.00209 ± 0.00045 

0.976687 

0.988829 



= 0.022442 


= 0.29 

= 0.31 

-0.00135 ± 0.00033 
-0.00112 ± 0.00034 

1.000000 

1.010363 

-0.00147 ± 0.00045 
-0.00116 ± 0.00046 

1.000000 

1.012431 


Table 2. Measured BAO peak positions qq and ot^ recovered from the mocks and presented in terms of their expected values. These results have been 
obtained from the QPM post-recon mocks for the LOWZ and CMASS samples, for different assumptions of cosmologies and different BAO models. We show 
the difference between the measured and the expected BAO position in the monopole, ao — Oq’'*’ and in the -moment, a 2 — Fof reference the 

expected BAO peak position for each cosmology and BAO model for both moments are listed. Two different cases are analysed: ct{P), where the BAO peak 
position is obtained from the fit to the mean of 1000 mocks; and a{P), where the BAO peak position is obtained from the mean of the BAO peak position fits 
to each individual mock. The former statistic is only sensitive to the systematic errors of the model, whereas the latter includes noise-dependent effects. 


by the potential differences between the true underlying Clm and 
the fiducial value assumed, in this case = 0.31. 

Fig. [^displays the output of the MCMC chains when estimat¬ 
ing the BAO peak position in the monopole and -moment, in the 
pre- and post-recon data catalogues, for the LOWZ and CMASS 
samples, as labeled in each panel. The blue and red contours display 
the 68% and 95.4% confident regions, respectively In this case the 
QPM covariance matrix has been used in all cases. By comparing 
the pre- and post-recon panels we see how the likelihood surface 
in the qq — a 2 parameter space is significantly reduced in the both 
LOWZ and CMASS sample, both for oq and a 2 - These results are 
consistent with DR 11 results jAnderson et al.|2014||Tojeiro et al.| 
|20T4l l. 

Table|^presents the marginalised measurements of qq and 0:2 
(and the derived AP parameters, ay and a±_ for the post-recon sam¬ 
ples) calculated from the MCMC chains (such as those presented 
in Fig.[^, for the following cases: 

(i) Pre-recon catalogues when the MD-Patchy covariance ma¬ 
trix has been used. 

(ii) Pre-recon catalogues adopting the QPM covariance matrix. 

(iii) Post-recon catalogues when the MD-Patchy covariance 
matrix has been used. 

(iv) Post-recon catalogues when the QPM covariance matrix has 
been used. 

The pre-recon results from both MD-Patchy and QPM co- 
variance matrices agree very well, both in the measurement and 
errors, suggesting that, as it was observed for the RSD analysis in 
Paper I, the impact of the differences between mocks when trans¬ 
lated into the covariance matrices of the power spectra is not sig¬ 
nificant. For pre-recon, we measure the BAO peak position in the 
monopole with ~ 3% and ~ 1.5% precision in the LOWZ and 
CMASS samples, respectively. In the -moment we are able to 
measure the BAO peak position with ~ 4% and ~ 2.4% precision 
in the LOWZ and CMASS samples, respectively. 

The post-recon results using either MD-Patchy or QPM co- 
variance matrices also agree very well, with differences that are 


Pre-recon 


Post-recon 



0.94 0.96 0.98 1 1.02 1.04 0.96 0.97 0.98 0.99 1 1.01 1.02 



Figure 2. Likelihood surfaces in the ciq — Q 2 parameter space obtained 
from measuring the power spectrum DR12 data monopole and ^^-moment. 
The black points are the (down-sampled) output of the MCMC chain, for 
the CMASS and LOWZ samples, for both pre- and post-reconstructed cata¬ 
logues, as labeled. In all cases the QPM-covariance matrix has been used for 
estimating the power spectrum errors and their correlations. The blue and 
red contours show the 68% and 95.4% confident regions, respectively. 


much smaller than the uncertainties due to the statistical errors. The 
post-recon results show an improvement on the level of precision 
for the BAO peak position determination, both in the monopole and 
in the -moment. We are able to improve the determination of the 
BAO peak position in the monopole down to ~ 1.8% and ~ 0.90% 
precision for the LOWZ and CMASS samples, respectively. The 
BAO peak position in the -moment is determined with ~ 2.8% 
and ~ 1.4% precision for the LOWZ and CAMSS samples, respec¬ 
tively. For the CMASS sample, this result represents an improve¬ 
ment of ~ 50%, both in ao and 0 : 2 . For the LOWZ sample the 
improvement is of ~ 40% for qq and ~ 25% for 02 . 
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Sample 

Catalogue 

Covariance 

ao 

0^2 

ro2 

“II 

ax 

r-yx 

X^/dof 

LOWZ 

Pre-recon 

MD-Patchy 

QPM 

0.998 ±0.028 
0.989 ±0.031 

1.019 ±0.037 
1.005 ±0.039 

0.78 

0.78 

1.053 ±0.067 
1.031 ± 0.069 

0.974 ± 0.039 
0.971 ± 0.042 

-0.33 

-0.32 

37.16/43 

35.72/43 


Post-recon 

MD-Patchy 

QPM 

1.009 ±0.017 
1.009 ±0.019 

1.018 ±0.027 
1.016 ±0.029 

0.81 

0.81 

1.032 ±0.050 
1.027 ± 0.053 

0.999 ± 0.023 
1.001 ± 0.025 

-0.40 

-0.42 

59.53/43 

55.57/43 

CMASS 

Pre-recon 

MD-Patchy 

QPM 

1.002 ±0.015 
0.997 ± 0.016 

0.996 ± 0.024 
0.988 ± 0.024 

0.75 

0.73 

0.987 ± 0.045 
0.974 ± 0.045 

1.010 ±0.023 
1.010 ±0.025 

-0.52 

-0.52 

43.73/38 

38.48/38 


Post-recon 

MD-Patchy 

QPM 

0.9895 ± 0.0091 
0.9899 ± 0.0088 

0.974 ± 0.013 
0.974 ± 0.014 

0.75 

0.74 

0.950 ± 0.024 
0.950 ± 0.025 

1.010 ±0.014 
1.010 ±0.014 

-0.43 

-0.51 

49.00/38 

37.24/38 


Table 3. Marginalised measurements of the BAO peak position for the monopole, and for the /i^-moment, cx. 2 , in the LOWZ and CMASS sample of the 
DR12 data as indicated. Also listed are their cross-correlation coefficient, ro 2 . For the post-recon samples we also show the derived AP parameters, a|| and 
a with their corresponding cross-correlation coefficient, r|| ^. We include the value of the best-fitting divided by the number of degrees of freedom (dof). 
The different rows are for the pre and post-recon catalogues, when both QPM and MD-Patchy covariances are used. 


The top panels of Fig. show the distribution of points in our 
MCMC chains for the post-recon LOWZ and CMASS DR12 data 
samples. We have included Icr (Ax^ = 2.30 in blue ellipses) and 
2(7 (Ax^ = 6.17 in red ellipses) contours from the Gaussian fit 
corresponding to the parameters extracted from the QPM covariance 
quoted in Table|^ For the BAO peak position variables, ao and Q 2 , 
as well as the AP variables ay and a±, the distribution of points 
matches well the Gaussian fits. 

Combining the best-fitting results on ay and ax of Table 
for both QPM and MD-PATCHY-based covariance matrices and the 
relations of Eq. we report, H{zi^owz)rs{zd) = (11.60 ± 
0.60) • 10^kms"'^and DA(«LOWz)/fs(2d) = 6.66 ± 0.16 for 
the LOWZ sample with a cross correlation coefficient of = 

0.41; and H{zcMASs)ra{zd) = (14.56 ± 0.37) • lO^kms”^ 
and DA{zcMASs)/'rs{zd) = 9.42 ± 0.13 with a cross corre¬ 
lation coefficient of thda ~ 0.47. These quantities represent 
for H{z)ra{zd) a 5.2% and a 2.5% measurement precision, for 
LOWZ and CMASS, respectively; and for DA{z)/ra{zd) a 2.4% 
and a 1.4% measurement precision, for LOWZ and CMASS, re¬ 
spectively. We can also put constraints on Dv{z) = {cz{l + 
z)D\H-^Y/^, which is the combination of parameters upon 
which the BAO location would depend for a galaxy sample with 
an isotropically distributed set of equally weighted pair separa¬ 
tions. We find that for LOWZ Dv( 2 L 0 wz)/»'s(zd) = 8.62±0.15; 
whereas for CMASS Dv{zcMASs)/rs{zd) ~ 13.70 ± 0.12. 


6.1 Consensus Values 


We combine the measurements presented above with those 
obtained using the correlation function monopo le, , and 


quadrupole, of the same dataset reported in iCuesta et al. 


|2016| companion paper), and report the consensus values for 
H{z)rs{zd), DA{z)/rs{zd^Sindi Dviz)/rs{zd). We summarise 
these values in Table|^ We see that both analyses agree well within 
1(7 error-bars. Although in theory the power spectrum and the cor¬ 
relation function contain the same information, in practice the anal¬ 
ysis is performed on a finite range of scales, and the information 
content of the two estimators differs slightly. The correlation be¬ 
tween them is large, as reported in the DR 11 analysis jAnderson| 
|et al.||20T4l> with a correlation factor of r ~ 0.95 for ao. Using 


^ Note that in jCuesta et al.|201^ companion paper) a different fiducial 
cosmology has been used to compute (ry and ax- However, the results 
in terms of H{z)re{zd) and DA{z)/rs(zd) should not depend on this 
assumption. 


the post-recon QPM mocks we have determined that for DR 12 data 
the correlation factor between the power spectrum and the correla¬ 
tion function is 0.91 in the LOWZ and 0.83 in the CMASS sample, 
for both ay and a±. This change is not produced by the differ¬ 
ences between the geometries of DRl 1 and DR 12, which are very 
similar, but because for DRll the fitting to the isotropic correla¬ 
tion function was performed using 8 bin-positions, whereas for the 
anisotropic correlation function of DR 12 was performed only using 
1 bin-position. This makes the DR12 anisotropic correlation func¬ 
tion fits noisier with respect to the DRl 1 ones, but at the same time, 
the correlation factor between ^ and P is reduced. In this way, the 
reduction on the error-bars produced by using more bin-positions is 
partially canceled by the fact that the correlation factor approaches 
to 1. Thus, the final results on the consensus value of the error-bars 
does not strongly depend on the number of bin-positions used in 
the correlation function analysis. 

From Table|4]we also note that both for LOWZ and CMASS 
the errors on H{z)rs{zd) are ~ 25% smaller when this quantity is 
estimated from the power spectrum. On the other hand, the errors 
on DA{z)/rs{zd) are similar for both LOWZ and CMASS. We 
have observed these differences not only from the likelihood of the 
data, but from the dispersion of the mocks, which demonstrates that 
the origin of this discrepancy is related to the methodology of how 
P™(fc) andP(^ ) are fit compared to how and are, and 
not with the data. These differences suggest that fitting the BAO 
peak position in the -moment (instead of fitting the quadrupole) 
is able to constrain better the AP distortion in the radial direction, 
both for LOWZ and CMASS. 

Using this correlation factor, we combine the measurements 
to obtain a consensus value. 


Qicon 


^ con 


Qp + Q;^ 

2 ’ 

(7p + CTJ /I + r \ 

J 


( 12 ) 

(13) 


where ap ± crp and a^ ± cr^ are the AP parameters measured 
from the power spectrum and correlation fraction, respectively, and 
aeon ± (7con the consensus value. When r = 0, both measurements 
would be independent and we would get a reduction of a factor 
of 2 on the error-bars, whereas when r = 1 both measurements are 
fully correlated and therefore the error-bars is unchanged. Using the 
correlation factors found in this paper, combining the correlation 
function and power-spectrum measurements improve the error bars 
by a few percent, both in LOWZ and CMASS. 
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Figure 3. Distribution of the MCMC points for the DR12 data post-recon catalogues where the QPM covariance matrix is used, in terms of the BAO peak 
position variables (oq ^ 2 , in the top panels) and the AP parameters («|| and a^, in the bottom panels) for the LOWZ sample (left panels) and for the 
CMASS sample (right panels). In blue lines the Icr ellipses (Ax^ = 2.30) and in red lines the 2cr ellipses (Ax^ = 6.17) corresponding to Gaussian fits to 
the likelihood based on the parameters given Table[^ For all cases the distribution of points is close to a Gaussian distribution. 


Sample 

Statistic 

[10® kms 

DA{z)/ra{zj) 

^hDa 

Dviz)/rs{zd) 


Power Spectrum 

11.60 ±0.60 

6.66 ±0.16 

0.41 

8.62 ± 0.15 

LOWZ 

Correlation Function 

11.65 ± 0.81 

6.67 ± 0.14 

0.29 

8.59 ±0.15 


Consensus 

11.63 ±0.69 

6.67 ±0.15 

0.35 

8.61 ± 0.15 


Power Spectrum 

14.56 ±0.37 

9.42 ± 0.13 

0.47 

13.70 ±0.12 

CMASS 

Correlation Function 

14.75 ± 0.50 

9.52 ± 0.13 

0.57 

13.79 ±0.14 


Consensus 

14.67 ±0.42 

9.47 ± 0.12 

0.52 

13.74 ±0.13 


Table 4. Distance constrain parameters, H{z)rs{z^), Da(AA and their correlation coefficient thDa’ inferred from the post-recon power spectrum 
(this work) and comelation function analysis jCuesta et al.|2016l companion paper), and the corresponding consensus value. The correlation coefficient between 
the power spectrum and the correlation function has been determined using the QPM post-recon mocks. In addition, we show the derived parameter Dy- 


1 CONCLUSIONS 


We have presented a new method to use Baryon Acoustic Os¬ 
cillations, measured with respect to the line-of-sight (LOS), to 
make cosmological measurements. The method fits the model of 


Anderson et al. 


j2014| to both the monopole and -moments. 


providing two correlated measurements of BAO positions using 
both statistics. This technique allows the clean separation of the 
BAO component from the broad-band power for both spherically- 
averaged determinations and measurements that are anisotropic 
around the LOS. This separation is harder to achieve when fitting 


the quadrupole where the BAO signal is removed if the fiducial 
cosmology matches the true one. We have applied this method to 
measurements of the anisotropic power spectra for the LOWZ and 
CMASS DR12 galaxies of the Baryon Oscillations Spectroscopic 
Survey of the Sloan Digital Sky Survey III, presented in Paper I, to 
constrain the BAO peak position in the monopole (ao) and in the 
/i^-moment (a 2 )- 

We have tested potential systematics of the BAO model using 
the post-recon QPM mocks catalogues, assuming different cosmol¬ 
ogy parameters and different BAO features models. We are able to 
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constrain the BAO peak position in the monopole and -moment 
better than 0.3% accuracy. These values are several times smaller 
than the Icr statistical errors in the LOWZ and CMASS samples. 
Furthermore, the scatter in these values appears random, and a sig¬ 
nificant component is likely due to residual noise. These tests ex¬ 
plicitly demonstrate that the method is independent of the fiducial 
cosmology assumed when calculating the power spectrum, or the 
BAO model to be fitted to the data. 

We have used the MD-Patchy and QPM mocks in the pre¬ 
reconstructed and post-reconstructed catalogues to estimate the co- 
variance matrices and hence the error-bars in the BAO peak po¬ 
sition parameters. We have found no significant differences in 
the results. From the post-reconstructed DR12 data power spec¬ 
trum monopole and ^^-moment we have measured: aol^LOwz) = 


1.009 ± 0.018 and 02 (^Lowz) = 1.017 ± 0.028 with ro 2 =0.81, 
where anowz = 0.32; and oo(zcmass) — 0.9897 ± 0.0090 
and a 2 (zcMASs) = 0.974 ± 0.014 with ro 2 = 0.75, where 
2CMASS = 0.57. We report these measurements in terms of cos¬ 
mological parameters: H{zi,owz)rs{zd) = (11.60 ± 0.60) • 
10^ kms“^ and DA(2LOwz)/ns(2d) = 6.66 ± 0.16 with a cross¬ 
correlation coefficient of thd^ = 0.41, for the LOWZ sam¬ 
ple; and JT(2cMASs))"s(2d) = (14.56 ± 0.37) • 10®kms“^ and 
TlA(zcMASs)/fs(-2d) = 9.42 ± 0.13 with a cross-correlation co¬ 
efficient of VHDji = 0.47, for the CMASS sample. 

We have combined these measurements with those obtained 
using the correlation function monop ole, and quadr upole. 


of the same dataset reported in (Cuesta et al. 


2016 


panion paper). We report H{zi,owz)rs{zd) = (11.63 ± 0.69) • 
10^ kms“^ and DA{zhowz)/rs{zd) = 6.67 ± 0.15 with a cross¬ 
correlation coefficient of thda = 0.35, for the LOWZ sam¬ 
ple; and -H'(zcMASs)f's(zd) = (14.67 ± 0.42) • lO^kms”^ and 
7lA(.2CMASs)/fs(:*ti) = 9.47 ± 0.12 with a cross-correlation co¬ 
efficient of thda ~ 0.52, for the CMASS sample. We see that the 
results reported from the analysis of the power spectrum agree well 
within 1(7 error-bars with the consensus values inferred from com¬ 
bining the power spectrum and 2-point correlation function analy- 
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